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ABSTRACT 

We use magnetohydrodynamical simulations of converging flows to investigate the 
process of molecular cloud formation and evolution out of the magnetised ISM. Here, 
we study whether the observed subcritical HI clouds can become supercritical and 
hence allow the formation of stars within them. To do so, we vary the turbulent Mach 
number of the flows, as well as the initial magnetic field strength. We show that 
dense cores are able to build up under all conditions, but that star formation in these 
cores is either heavily delayed or completely suppressed if the initial field strength 
is B > 3 fiG. To probe the effect of magnetic diffusion, we introduce a tilting angle 
ip between the flows and the uniform background magnetic field, which mimics non¬ 
ideal MHD effects. Even with highly diffusive flows, the formed cores are devoid of 
star formation, because no magnetically supercritical regions are build up. Hence we 
conclude, that the problem of how supercritical cloud cores are generated still persists. 

Key words: magnetohydrodynamics (MHD) - turbulence - ISM: clouds - ISM: 
kinematics and dynamics - ISM: magnetic fields 


1 INTRODUCTION 

Stars and stellar systems form within the densest regions of 
molecular clouds, in gravitationally unstable cores which re¬ 


side at the junctions of filaments (e.g. Andre et al. 2013 
|2014 ). Prior to gravitational collapse the build-up of fil¬ 


aments and the respective substructures is primarily con¬ 
trolled by magnetic fields and supersonic turbulence (e.g., 
Shu et aT7||198T[ |Mac Low &; Klessen||2004| |Crutcher et al.| 


2010| . But, the importance of magnetic fields for star for¬ 
mation is still debated (see e.g., Li et al.|2 014; Padoan et al. 


2014). On the one hand, the idea of supersonic turbulence 


controlling the star formation process assumes less impor¬ 
tant magnetic fields and thus primarily supercritical states. 
In such a scenario, the magnetic field lines are dragged along 
with the flow and density enhancements will collapse as soon 
as they become Jeans unstable. Furthermore, the turbu¬ 
lence is then not only supersonic but also superalfvenic (e.g., 
Padoan et al.|1999 Padoan &; Nordlund|1999 ). This leads to 
highly twisted field lines and the resulting molecular clouds 
and clumps will not be coherent entities. The morphology 
instead will be influenced by the statistics/nature of the tur¬ 
bulence. 

On the other hand, Mestel &; Spitzer] ( |1956| first quantified 
the influence of magnetic fields on star formation by intro- 
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ducing the mass-to-magnetic flux ratio fi = M /$ as a mea¬ 
sure of the relative importance of gravitational and magnetic 
energies. Usually, this quantity is normalised to its critical 
value ii c ~ 0.13 /y/G (or ii c ~ 0.1 6/VG for more sheet-like 
clouds(Nakano & Nakamura ) 1978| ). If the magnetic field is 
strong enough, accretion onto the cloud complex is medi¬ 
ated by the Lorentz force and mainly parallel to the field 
lines (e.g. Kudoh et al.|2007 Inoue &; Inutsuka|2008 Kudoh 


& Basu 2010; Hennebelle 2013). In the cloud interior, strong 


fields stabilise the filaments and clumps against gravity. This 
also results in a reduced fragmentation efficiency. 
Observationally, it has been shown in recent years that the 
magnetic field is indeed crucial for the star formation pro- 


cess (e.g., 

Beck 

2001 Crutcher et al. 

2010 

Li et al. 

2010 

Crutcher 2012 

Li et al. 2014 Pillai et al. 

2014). Li et al. 

(2010) used sub 

-mm polarisation measurements to retrieve 


the morphology of the magnetic field in molecular clouds 
and Galactic spiral arms. The authors have shown that the 
overall morphology of the field does not change significantly 
from the large scales down to the inner parts of molecular 
clouds. By using HI, OH, and CN Zeeman measurements 


Crutcher et al. (2010 see also Crutcher (2012),fig. 7) have 
shown that nearby molecular clouds and cloud cores can 
be separated into two regimes according to their column 
density and magnetic field strength. At low column densi¬ 
ties, the magnitude of the field almost does not change for 
roughly two orders of magnitude with a median value of 

















































































2 Kortgen & Banerjee 


Blos ~ 5 — 6 fiG. This regime also coincides with magneti¬ 
cally subcritical HI clouds. As was pointed out by [Crutcher] 
(2012 and references therein), these data are primarily dif¬ 
fuse HI clouds that are not self-gravitating , but are rather 
in pressure equilibrium with their surroundings. Note, that 
the total magnetic field strength will even be larger. E.g. 
recent studies give average values of 5 — 15/iG (Beck 2001 
Cr utcher et al.| 2010). 


At higher column densities the field strength increases close 
to linear, with increasing column density. At this stage, al¬ 
most all measurements indicate (super-)criticality. Conver¬ 
sion of these data points to volume density shows that the 
scaling in the latter region is B oc n 0 ' 65 , which perfectly fits 
to conditions of frozen-in magnetic field lines and isotropic 
collapse (which gives a relation B oc n 0-66 ). 

Numerically, the issue of magnetic fields and their relevance 
for molecular cloud formation has been investigated by many 
authors (e.g., Inoue V Inutsuka |2008| [Price &; Bate||20d8 


2009| |Kudoh V Basu||2010| | Vazquez-Semadeni et al.||2011 


Inoue &; Inutsuka 2012 Chen V Ostriker 2014). Most of 


them concentrated on the initial stages of the formation pro¬ 
cess. Already at this early temporal stage, the magnetic field 


was shown to be crucial. Price V Bate (2008) conducted 


smoothed particle hydrodynamics simulations of a 50 Mq 
molecular cloud of radius R = 0.375pc including magnetic 
fields of different strength (parameterised by critical mass- 
to-flux ratios of oo, 20, 10, 5, and 3). They found that strong 
fields tend to suppress fragmentation on the one hand and 
the formation of stars on the other hand. However, as they 
point out, strong fields generate voids within the molecular 
cloud, which are magnetically supported with plasma-/^ > 1. 
In addition, on very small scales, magnetic tension is able 
to prevent multiple fragments from merging, thus promoting 
fragmentation. 

More consistent with this study is the work by |Heitsch| 
et al. (2009). They have used MHD simulations of converging 


flows to analyse the impact of magnetic field strength and 
orientation on the formation of (molecular) clouds. Specifi¬ 
cally, they looked at the extreme cases of the magnetic field 
being either aligned with or perpendicular to the flow direc¬ 
tion. The flows were driven continuously due to the choice 
of inflow boundary conditions. Hence, the mass-to-flux ra¬ 
tio in their study would approach infinity in the limit of 
infinite timescales. Note that the authors have not included 
self-gravity in their simulations. Thus, every overdense sub¬ 
structure is pressure confined. However, they identify fila¬ 
ments and clumps that form due to turbulent compression, 
with clouds becoming more filamentary if magnetic fields are 
included. But, it is the alignment of the magnetic and (ini¬ 
tial) velocity field that controls the formation of dense struc¬ 
tures. As the authors point out, clouds are able to condense 
out of the WNM, if the fields are aligned. In case the mag¬ 
netic field is perpendicular to the inflows, magnetic pressure 
suppresses the formation of dense structures, which could 
be termed molecular. However, there exist regions, which 
merge to form a filamentary network of diffuse gas. 


Inoue V Inutsuka (2009) studied the evolution of the shocked 


slab between two converging flows in the ISM by means of 
two-fluid MHD simulations in a 30 pc x 10 pc box. The au¬ 
thors varied the angle between the mean magnetic field and 
and the flows. From analytical estimates they found a crit¬ 
ical velocity, which depends on the magnetic field strength 


and the mentioned angle. If the flow velocity is larger than 
the critical velocity only HI clouds are able to form because 
of dominating magnetic pressure. If it is less than the critical 
velocity, dense molecular clouds condense out of the WNM 
within the shocked slab. As the authors also point out, the 
dependence on the angle is crucial for the evolution of the 
gas within the slab, since the critical velocity goes to zero 
for angles approaching 90°. 


Most recently, Chen V Ostriker (2014) studied the formation 


of prestellar cores due to the convergence of gas flows within 
molecular clouds. In detail, their simulation box was about 
1 pc, representing a collapsing molecular clump. In order to 
analyse the core formation process, they used ideal MHD 
as well as non-ideal MHD via ambipolar diffusion (AD). In 
all of their models core formation was initiated by the col¬ 
lision of gas streams along the background magnetic field. 
With AD only the later stages were seen to differ from the 
ideal MHD models, since the density regimes where AD is 
becoming efficient are build up via accumulation of gas by 
colliding flows. The mass-to-flux ratio of the cores formed 
in their simulations is in the range /i/ti c ~ 0.5 — 7.5 with 
a median value of fi/pL c ~ 3. Thus, most of the cores are 
supercritical 

However, it is important to conduct large scale simulations 
in order to take into account the whole evolutionary track 
of the gas from the diffuse ISM to the dense cores. This was 


achieved by Vazquez-Semadeni et al. (2011), who analysed 


molecular cloud formation subject to magnetic fields of dif¬ 
ferent initial strength. It has been shown that stronger fields 
tend to delay the onset of star formation. The authors (see 
also Hartmann et al. 2001; Vazquez-Semadeni et al. 2006, 


2007) point out that the diffuse gas becomes molecular, self- 


gravitating and magnetically supercritical at the same time. 
This simple approach can explain the subcriticality of the 


diffuse HI clouds shown in Crutcher et al. (2010); Crutcher 


(2012). Heitsch V Hartmann (2014) mention that the su¬ 


percritical state can be reached via gas accretion along the 
magnetic held lines. But as was stated by |Hartmann et ah] 
(2001), the accumulation length to become magnetically su¬ 


percritical is 


! 470 


( Bo ) 

( U i 

V nGJ 

V 1 cm -3 ) 


pc. 


a) 


Vazquez-Semadeni et al. (2011) argue that, since the mag¬ 


netic held lines in the Galactic plane describe closed circles, 
this length scale is easily overcome. This also indicates that 
the mass-to-hux ratios are lower limits and the data points 


shown in Crutcher et al. (2010) are only a temporal stage of 


subcriticality. But, as Carroll-Nellenback et al. (2014) point 


out, how lengths of L > 500 pc are too large in order to sus¬ 
tain a large scale coherent how. Bulk motions of this order of 
magnitude should rather fragment due to supersonic turbu¬ 
lence and thus diminish. Hence, the build-up of supercritical 
clouds would be delayed or even suppressed completely. The 
process, how molecular clouds achieve the transition from 
sub- to supercritical states is thus still an open question. 


1 Note that the authors use inflow boundary conditions for the 
two converging hows. Hence, the mass-to-hux ratio of the simu¬ 
lation domain will grow with time and so it will for the cores as 
they accrete mass from an practically inhnite mass reservoir. 

























































































MCs and magnetic fields 3 


In this study we tie in with the work of |Vazquez-Semadeni 


et al. (2011) by determining molecular cloud formation un¬ 


der different initial conditions. In section [2] we therefore 
introduce our numerical model and the initial conditions. 
Section |3] deals with the formation and evolution of clouds 
formed by head-on converging WNM streams under varying 
initial conditions. The following section [4] then introduces 
the tilt of one flow with respect to the magnetic held and 
discusses in detail the evolution of the clouds and their sub¬ 
sequent star formation activity. This study is closed by a 
brief summary in section [5] 


2 NUMERICAL SETUP AND INITIAL 
CONDITIONS 

2.1 Details of the numerics 

For this study we use the finite volume AMR code FLASH 
(Fryxell et al. 2000| Dubey et al. 2008). During each timestep 
a Riemann problem is solved at the cell interfaces, yielding 
the respective fluxes for the hyperbolic partial differential 
equations. The MHD fluxes are computed by a multi wave 


Riemann solver developed by Bouchut et al. (2007 2009) 
and implemented in FLASH by Waagan et al. f( |2011| ), which 
preserves positive states for density and internal energy. 

We apply periodic boundary conditions for the (magneto- 
hydrodynamics and isolated ones for gravity. In addition to 
the basic (ideal) MHD equations, we include selfgravity as 
well as heating and cooling. The latter is treated as a source 
term in the energy equation and we follow the recipe by 
Koyama V Inutsuka (2000, with modifications by Vazquez- 
Semadeni et al. (2007)) for the radiative cooling of the gas. 
Since we are interested in the process of star formation, we 
also include sink particles to follow truely collapsing regions 
(Fed errath et al.|2010 ) and the local Jeans length is resolved 
with at least ten grid cells to fulfill the Truelove criterion 
( Truelove et al.|1997 ). In order to replace a certain gas vol¬ 
ume by a sink particle, the gas has to pass several checks, 


which are described in great detail in Federrath et al. (2010). 


We use a density threshold of n ~ 2 x 10 b cm - A Once a sink 
particle has been created, it is only allowed to accrete gas. 
Feedback is not included. 

The numerical grid is refined when the local Jeans length 
is resolved with less than ten grid cells and derefmed if it 
consists of more than 100 cells. 


2.1.1 Treatment of Ambipolar Diffusion 

We have conducted one simulation including the non-ideal 
MHD effect of ambipolar diffusion (AD). The AD module 
was implemented in FLASH and extensively tested by |Duf-| 
fin V Pudritz (2008). It uses the strong coupling approxima¬ 
tion (like |Chen fc Ostri ker (2012, see also Chen V Ostriker 


( 2014| )), which was shown to be valid in the physical regime 


we are analysing (see appendix in Vazquez-Semadeni et al. 
2011| . However, since we are using a slightly different den¬ 


sity threshold than in Vazquez-Semadeni et al. (2011) and 


different turbulent Mach numbers, the validity has to be 
proven again: 

Taking a typical length scale of l = 0.0625 pc (which corre¬ 
sponds to the accretion/softening radius of the sink particles 



Figure 1. Setup of the initial conditions. The black dashed lines 
denote the flow axes and the position of the collision layer, re¬ 
spectively. Adapted from [Vazquez-Semadeni et al~| ( |2007| . 


in our simulation with 11 levels of refinement) and a typical 
velocity at these scales of 0.5 km/s (taken from fig. [ 5 ]) the 
ratio 


mI/Rad(1 ) « 1.5 x 10" 


( 2 ) 


Here Ma is the Alfven Mach number and Rad( 0 is the AD 
Reynolds number at scale l. Hence, according to |Li et al.] 
(2006, see also Vazquez-Semadeni et al. ( 2011) )) the strong 
coupling approximation is satisfied very well in our simula¬ 
tions. 

The fluxes are computed using a central differencing scheme 
and the numerical timestep is primarily controlled by AD. 


Note that the implementation by Chen V Ostriker (2012) 
uses super-timestepping to speed up the simulation (see also 
Choi et al. 2009). 


2.2 Initial conditions 


Our numerical setup is very si milar to this of [V azquez- 


Semadeni et al. (2007 see also Vazquez-Semadeni et al. 


(2011) and Banerjee et al. (2009)). The physical size of the 
numerical box is 256 pc and the two cylindrical flows col¬ 
lide at the centre of the domain, that is, at L x = 128 pc. 
Each flow has a linear dimension of l = 112 pc and a ra¬ 
dius of r — 64 pc, thus twice as large as in the aforemen¬ 
tioned studies (see table [l] for the runtime details). The 
heating and cooling prescription by |Koyama h Inutsuka| 
(2000) gives a thermally unstable regime in the density range 
lcm -3 ^ n ^ 10 cm -3 , corresponding to a temperature in¬ 
terval of 500 K ^ T ^ 5000 K if thermal equilibrium con¬ 
ditions are applied. According to these restrictions we have 
chosen the initial density to be n = 1 cm -3 and the temper¬ 
ature as T — 5000 K. With this temperature, we can define a 
sound speed of the warm neutral medium and we choose the 
velocity of every single flow in such a way that the result¬ 
ing isothermal sonic Mach number is Atf = 2. In addition 
we add a turbulent velocity held to the hows to mimic the 


general turbulent behavior of the ISM (Mac Low & Klessen 


256 pc 













































































































4 Kortgen & Banerjee 


Table 1. Overview of the conducted simulations with varying flow and magnetic field parameters. (f> is the inclination of one of the flows 
and A4f is the isothermal Mach number of the converging WNM streams. Mrms denotes the Mach number of the turbulent fluctuations, 
A4a is the turbulent Alfven Mach number, f} is the ratio of thermal to magnetic pressure and k* indicates the ratio of turbulent to 
numerical diffusion with the turbulent diffusion coefficient calculated according to |Lazarian et al. (2012) (their equations (3) and (4)). 
Numerical diffusion is evaluated as ft n um = Ax x va x , where v/±x is the velocity at grid scale Ax. % denotes the effective ratio of magnetic 
to numerical diffusion (see appendix 0- /jl//i crit is the normalised mass—to-magnetic flux ratio. 


Run Name 

<t> 

(°) 

in 

(/xG) 


Mrms 

Ma 

K* 

X 

p 

M/Mcrit 

(pc) 

Min. Ax 6 ) 

B3M0.4I0 

0 

3 

2 

0.4 

0.39 

2.53 

0.00 

1.93 

0.79 

0.03 

B3M0.8I0 

0 

3 

2 

0.8 

0.79 

21.04 

0.00 

1.93 

0.79 

0.03 

B3M1.2I0 

0 

3 

2 

1.2 

1.18 

42.67 

0.00 

1.93 

0.79 

0.03 

B4M0.4I0 

0 

4 

2 

0.4 

0.29 

1.04 

0.00 

1.08 

0.59 

0.03 

B4M1.5I0 

0 

4 

2 

1.5 

1.10 

42.67 

0.00 

1.08 

0.59 

0.03 

B5M0.5I0 

0 

5 

2 

0.5 

0.29 

1.04 

0.00 

0.69 

0.47 

0.0075 

B3M0.5I30 

30 

3 

2 

0.5 

0.49 

5.02 

6.91 

1.93 

0.79 

0.0075 

B3M0.5I50 

50 

3 

2 

0.5 

0.49 

5.02 

16.22 

1.93 

0.79 

0.0075 

B3M0.5150a 

50 c ) 

3 

2 

0.5 

0.49 

5.02 

16.22 

1.93 

0.79 

0.0075 

B3M0.5I60 

60 

3 

2 

0.5 

0.49 

5.02 

20.73 

1.93 

0.79 

0.0075 

B3M0.8I60 

60 

3 

2 

0.8 

0.49 

5.02 

20.73 

1.93 

0.79 

0.0075 

B4M0.5I30 

30 

4 

2 

0.5 

0.36 

1.99 

9.22 

1.08 

0.59 

0.0075 

B4M0.5I60 

60 

4 

2 

0.5 

0.36 

1.99 

27.63 

1.08 

0.59 

0.0075 

B5M0.5I30 

30 

5 

2 

0.5 

0.29 

1.04 

11.52 

0.69 

0.47 

0.0075 

B5M0.8I30 

30 

5 

2 

0.8 

0.47 

4.43 

11.52 

0.69 

0.47 

0.0075 

B5M0.5I40 

40 

5 

2 

0.5 

0.29 

1.04 

19.03 

0.69 

0.47 

0.0075 

B5M0.5I50 

50 

5 

2 

0.5 

0.29 

1.04 

27.03 

0.69 

0.47 

0.0075 

B5M0.5I60 

60 

5 

2 

0.5 

0.29 

1.04 

34.55 

0.69 

0.47 

0.0075 

B5M0.8I60 

60 

5 

2 

0.8 

0.47 

4.43 

34.55 

0.69 

0.47 

0.0075 

B5M0.5I60Mf4 

60 

5 

4 

0.5 

0.29 

1.04 

34.55 

0.69 

0.47 

0.0075 

B5M0.5I60AD d ) 

60 

5 

2 

0.5 

0.29 

1.04 

34.55 

0.69 

0.47 

0.0075 


Remarks: 

а) According to the prescription by 

б) Maximum allowed resolution i: 


Nakano & Nakamura) |l978| (i.e. p cr ^ ~ 0.1 6/y/G). 
re simulations. 


c) Simulation with a different initial turbulent seed field. 

d ) Run with ambipolar diffusion. Simulation was stopped at t « 12 Myr. 


2004). The turbulent fluctuations are calculated in Fourier 
space with a Burger’s type spectrum, i.e. E(k) oc k~ 2 for 
k ^ fonj, where fo n j is the energy injection scale. Further¬ 
more, these turbulent fluctuations trigger the onset of dy¬ 
namical instabilities such as the non-linear thin-shell insta¬ 
bility (NTSI, (|Vishniac||1994|) or the Kelvin-Helmholtz in¬ 
stability (e.g. Heitsch et al.||2005 2008). The initially uni¬ 
form magnetic field has a strength of \B] = {3,4, 5} /jlG and 
is aligned with the flows, that is, B oc x where x is the 
unit vector in the x direction. Using the description for the 


critical mass-to-flux ratio by Nakano V Nakamura (1978), 
i.e. /icrit = 0.16/v^S, the two streams in total are initially 
subcritical, but can become supercritical very fast due to 
accretion of mass along the field lines (see table [l]) . The 
numerical resolution is adjusted to give a maximum refine¬ 
ment level of £max = 13, which corresponds to a maximum 
physical resolution of Ax max ~ 0.007 pc « 1600 AU. 


MC FORMATION BY HEAD-ON 
COLLIDING FLOWS 


Colliding streams of gas are ubiquitous in the ISM (e.g. due 
colliding supernovae shells, Inoue & Inutsuka 2008) as well 
as in the interior of molecular clouds (e.g. in filaments or 
the junctions of filaments, Hennebel le et al.||2008| |Chen 
Ostriker 2014). Therefore the dynamics and the structure 


can vary significantly, depending on the galactic or local 


environment or the respective driving mechanism (e.g. In¬ 
oue & Inutsuka 2008, 2012). In this section we summarise 


the evolution of molecular clouds, which are being formed 
by head-on colliding flows. For more thorough analyses, 


(2009); 

Vazquez-Semadeni et al. 

(2006); 

Vazquez-Semadeni 

et al. (2011); Hennebelle &; Perault (1999); Hennebelle et al. 


physical parameters of the respective simulations is given in 
table [l] 



































































MCs and magnetic fields 5 




Figure 2. Left: Evolution of cloud (solid) and sink particle (dashed) masses for three turbulent velocity fields (for a colour version see 
online manuscript). The onset of star formation is clearly seen to be delayed due to the action of turbulent motions that keep dissolving 
dense structures. The horizontal dash-dotted, black lines denote the masses of observed molecular clouds, like the Perseus MC (PMC, 
|Lombardi et al.|2010| ), the Taurus MC (TMC, [Lombardi et al.|2010| ), and the Corona Australis complex (CA, [Alves et al.|2014| ).Middle: 
Star formation efficiency for B3M... I... runs. Right: Corresponding root mean square velocity of the dense gas. The RMS-velocities 
converge after turbulence has decayed and global collapse of the cloud has begun. Before this point, the amplitude of the resulting 
turbulent velocities is determined by the initial conditions. 
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Figure 3. Slices of column density along two different axes for different simulations at t^lOMyr. Top: Along x-axis (i.e. face-on view). 
Bottom: Along y-axis (i.e. edge-on view). 


3.1 Varying the turbulent velocity 

The dynamics of molecular clouds which formed in the com¬ 
pression zone of two colliding streams strongly depend on 
the initial kinematics of the individual flows. On the one 
hand, the flows are supersonic with respect to the WNM 


and thus generate strong shocks and compressions (Vazquez- 


Semadeni et al. 2007 Banerjee et al. 2009). On the other 


hand, large scale instabilities as well as stellar feedback in¬ 
ject energy into the ambient ISM. This energy, if not al¬ 
ready in the form of kinetic energy, can be converted to 
kinetic energy and thus a turbulent regime is produced, 
where the turbulence cascades down until it is dissipated on 


atomic/molecular scales. This turbulence is primarily su¬ 
personic (e.g., Mac Low & Klessen 2004). These random 


motions generate a certain level of anisotropy within the 
bulk flows and the respective contribution to the process 
of molecular cloud formation is two-folded. Firstly, turbu¬ 
lence contributes an effective ram pressure, which can help 
to stronger compress fluid elements. Secondly, the inhomoge¬ 
neous velocity field distorts the overall bulk flow and reduces 
the mass flux, which then directly translates to the build up 
of less massive clouds (see fig. [2]) . 

If the collision of the WNM streams is along the magnetic 
field lines, the early stages (t ^ 3 Myr) of cloud formation 





















































































6 Kortgen & Banerjee 


can be understood as being nearly independent of the mag¬ 
netic field. The first phases during the collision are thus 
controlled by the bulk and turbulent velocity (see tab. [I]). 
Fig. [ 2 ] shows the evolution of the dense gas (n ^ 100 cm -3 ) 
for different initial turbulent Mach numbers. The compres¬ 
sion by the flows induces the formation of a molecular cloud 
by the combined action of dynamical and thermal instabil¬ 
ity (e.g., Vazquez-Semadeni et al.|2007 Heitsch et al. 


Banerjee et al.||2009 ). Due to the onset of runaway cooling 

of thermally unstable gas, the cloud becomes more massive 
with time. At the same time it assembles mass by accretion 
of gas along the field lines. Independent of the degree of tur¬ 
bulence, the onset of dense gas formation starts at the same 
time indicating the dominance of the ram pressure by the 
bulk flows. Only at slightly later times around f « 2 — 3 Myr 
the effects of different turbulent Mach numbers are seen. 
Flows of higher turbulent Mach numbers reduce the mass 
flux and hence reduce the final mass of the cloud. This is 
seen in fig. [3] The stronger the turbulence, the less compact 
is the resulting cloud. The mass concentrates in pressure 
confined filaments, which are further immersed in a diffuse, 
warm medium, with a steep density and temperature gradi¬ 
ent between the WNM and CNM that can be interpreted as 
a phase-transition front rather than a contact discontinu¬ 
ity due to the ambient mass flux across the transition layer 
( Banerjee et al.||2009 ). 

At later stages, the initial turbulence has decayed and the 
presence of turbulent motions is due to self-gravity (see 
e.g. Ballesteros-Paredes et al. 2007). Once self-gravity dom¬ 
inates, certain regions then proceed to collapse to form a 
star. The onset of star formation is clearly delayed by the 
presence of stronger initial turbulence (see fig. [2|. 


3.2 Dependence on the magnetic field strength 


In the previous section we have neglected the possible in¬ 
fluence of the magnetic field. However, the ISM and molec¬ 
ular clouds are highly magnetised (in the range 5 — 15/xG, 


Beck 

2001 1 

2010, 

2014) 


2014) and thus the magnetic field affects the overall 


evolution of molecular clouds in the ISM as well as their 


preceding condensation out of the latter (e.g., Hennebelle fe 
Perault|1999 Hennebelle|2013 ). As was shown by Crutcher 


et al. (2010) using Zeeman measurements, the line-of-sight 


component of the interstellar magnetic field can be approx¬ 
imated by an interval of nearly constant magnitude fol¬ 
lowed by a regime that consists of a linear increase of the 
field strength as function of (column-)density. Since Zeeman 
splitting provides information of one component only, the 
total magnetic field strength will be larger. Recent studies 
provide average values of the magnetic field of approximately 
5 — 15/xG ( Beck||2001 Crutcher et al.||2010 ). It is thus rea¬ 
sonable to investigate the influence of varying magnetic field 
strength on the molecular cloud formation process. This has 
recently been done by Vazquez-Semadeni et al.| (2011) for 
initial magnetic field strengths of \B\ — {2,3,4}/xG (corre¬ 
sponding to n = {1.18,0.79, 0.59}fj and the action of am- 


2 Note, the values for the mass-to-flux ratio in | Vazquez^] 
|Semadeni et al.| ( |2011| ) refer to the box length of 256 pc, instead 
of the flow length of 112 pc, which we here take care of. 



Figure 4. Evolution of cloud (solid) and sink particle (dashed) 
masses for three magnetic field strengths (for a colour version 
see online manuscript). The clouds assemble similar final masses, 
but some slight differences in the accretion phase are seen. More 
prominent is the lack/delay of star formation for simulations with 
\B] > 3 /xG. The horizontal dash-dotted, black lines denote the 
masses of observed molecular clouds, like the Perseus MC (PMC, 
Lombardi et al. ||2010| ) , the Taurus MC (TM C, |Lombardi et al.| 
2010J), and the Corona Australis complex (CA, | Alves et al.|2014^ 


bipolar diffusion. Our study covers the upper range of their 
values, namely the range of \B\ = {3,4, 5} /xG. The choice of 
these field strengths gives thermally dominated ( \B\ =3 /xG) 
environments, regimes with an equipartition of thermal and 
magnetic energies ( \B\ — 4/xG), and completely magnet¬ 
ically dominated regions ( \B\ = 5/xG). The evolution of 
the cloud and sink particle mass for different initial field 
strengths is shown in fig. Here, the final cloud masses do 
not differ too much from each other, showing that the ini¬ 
tial turbulent motions are more efficient in controlling the 
early phases of gas accumulation. But differences are seen 
in the early mass accretion. The cloud, which is embedded 
in a strong magnetic field is seen to be build up at a slightly 
earlier time after the start of the simulation. Furthermore 
the accretion of matter from the diffuse halo surrounding 
the cloud at early times differs for the strongest initial mag¬ 
netic field. This fact can be explained by momentum and 
energy conservation. Since the flows collide head-on, the gas 
is compressed in the collision layer. The external ram pres¬ 
sure by the flows forces the gas to move perpendicular to the 
magnetic field lines in order to ensure conservation of lin¬ 
ear momentum. At the same time gas compression enhances 
the magnetic field strength. Magnetic tension then acts as a 
restoring force and since the plasma -[3 is less than unity, the 
dominant magnetic field is too stiff to be bend efficiently. 
This results in a less efficient gas motion perpendicular to 
the original bulk flow motion and an earlier compression 
of the gas. Thus, efficient accretion happens only along the 
field lines. The rightmost plot in fig. [ 3 ] shows the column 
density after t ~ 10 Myr for a strong magnetic field. The 
density gradients are smoother in comparison to the weaker 
field and the cloud is more compact, i.e. no clearly defined 
filaments condense out. At the same time the column den- 
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sity (in the face-on view) does not reach sufficiently large 
values. The critical column density to become magnetically 


supercritical can be written as (see also Vazquez-Semadeni 
et al.|2011 1 


iVcrit « 2.92 x 1<T 


(±-V 

Wg 


(3) 


Although all different clouds assemble mass by accretion 
from the surrounding diffuse gas, the most striking differ¬ 
ence is the complete lack of star formation for the higher 
magnetised clouds. In order to quantify the gas dynamics 
at a specific evolutionary stage, fig. [5] shows the probability 
distribution function of the column density (left, hereafter 
N-PDF) and the mass-to-flux ratio (right, from now on /jl- 
PDF). As shown by Vazquez-Sema deni V Passot] pQQQ), the 
statistics of a gas can be analysed by using a density PDF. 
For isothermal turbulence, this PDF develops a lognormal 
distribution with its variance depending on the Mach num¬ 
ber of the gas. More recently it has been demonstrated that 
the width of the distribution also depends on the plasma- 
P as well as on the turbulent forcing parameter, i.e. if the 
driving of turbulence is purely solenoidal or compressive (see 
e.g., Federrath V Klessen 2012 2013). The same is also valid 
for the N-PDF, which is used in observational studies since 


the volume density is not accessible (e.g., Kainula inen et al. 
2011 Schneider et al.|2013 2014). The shape of the N-PDFs 


m fig. [5] is not lognormal. The reason is the multi-phase na¬ 
ture of the ISM. However, in the weaker magnetised case, a 
power-law tail at high column densities evolves, which is al¬ 


ways seen in self-gravitating systems (Federrath V Klessen 
2013 Schneider et al. 2013), indicating the presence of grav¬ 


itationally unstable regions. The vertical lines denote the 
threshold column densities according to eq. [ 3 ] The thermally 
dominated case shows a transition to supercritical states, 
whereas the maximum column density in the magnetically 
dominated gas is approximately a factor of five lower. At this 
time the WNM flows vanished. Increases in column density 
are only due to mass accretion from the environment. Since 
the column density and the mass-to-flux ratio are coupled 
via 


E 

| -Hlos I ’ 

the column density in g/cm 2 
, the overall shape of the /i-PDF should be very 
similar to the one of the N-PDF. This is indeed the case, 
as can be seen from fig. [5] right. The modifications are due 
to the additional dependence on the magnetic field. Here, 
again, the mass-to-flux ratio shows a similar distribution, 
but for the weaker field it is shifted towards higher values. 
Furthermore, the transition from subcritical to trans-critical 
regions is smoother, because of the lack of stiffness of the 
magnetic field. The dependence on column density then also 
implies the outcome of a power-law tail in the distribution, 
which continues up to values of log (fi/fi c ) ~ 1, showing the 
presence of highly unstable, dynamically dominated regions. 
The mass-to-flux ratio for runs B4M0.4I0 and B5M0.5I0 
is similar distributed. This indicates that initial dynami¬ 
cal processes should be more energetic than observed in 
the simulations, since dynamic compressions always result 
in increasing magnetic energy, which at some stage starts to 
dominate over thermal and gravitational energy. This yields 


where E is 
mura 1978b 


(4) 

(Nakano V Naka- 



Figure 6. Setup of the initial conditions with inclined flow. The 
box and physical parameters are the same as in fig. ^ but now 
one flow is inclined at an angle ip with respect to the background 
magnetic field. The inclined axis is shown as the red dashed line. 
Note, the collision still occurs at the centre of the simulation 
domain, in contrast to what is schematically shown here. 


a re-expansion of compressed regions and a simultaneous 
stabilisation of these. 

The results from simulations with higher magnetisation now 
raise the question, how stars can form in such highly mag¬ 
netised media. 


4 INCLINED WNM FLOWS 


Here we probe the influence of inclined colliding flows. In¬ 
clined collisions are easily justified by assuming the emer¬ 
gence of a supernova shock wave and its propagation through 
a Galactic spiral arm or by non-uniform large scale gravita¬ 
tional forces. The motion of the flow at an inclination with 
respect to the magnetic field results in an enhanced diffu- 
sivity of the latter and this process thus can be thought of 
as a non-ideal MHD process (e.g. Heitsch et al. 2005 Inoue 


& Inutsuka 2008; Heitsch et al. 2009) 


4.1 The setup 

The initial geometry can be seen in fig. [6] The basics are the 
same as for the head-on case, but now one flow is inclined at 
an angle ip with respect to the x-axis. The initial background 
magnetic field is kept constant and aligned with the x-axis. 
The figure may imply that there might be a region where 
quiescent gas resides, but this is not the fact. The two flows 
still collide in the centre of the simulation box and since 
the magnetic field is still uniform, the collision will induce 
a normal shock. Table [l] lists the initial parameters for this 
study. 


256 pc 
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iog(n/n c ) 

Figure 5. Left: Column density PDF for different initial magnetisations at 20Myr. The weaker magnetic field allows for the development 
of a power-law tail, which indicates gravitational collapse. The vertical lines denote the respective threshold column densities for magnetic 
criticality after eq. [3] Right: Corresponding /x-PDF. Here, the vertical line indicates the critical ratio. 


4.2 Magnetic Flux Reduction and Star Formation 

We commence with the thermally dominated case, i.e. 
\B_ 0 \ = 3 fiG. The inclination of one flow increases the diffu- 
sivity (see tab. [l] and appendix [Af where larger angles lead 
to larger magnetic diffusivity. Fig. [2] has already shown the 
evolution of the cloud mass as function of time. The mass 
crucially depends on the strength of the initial turbulent 
velocity fluctuations. In addition, if one applies an inclined 
flow, shearing motions and magnetic effects have to be taken 
into account. The magnetic field is able to slow down the 
inclined flow so that the collision will end soon and no gas 
is driven into the thermally unstable regime. But the final 
masses of the formed clouds are very similar, only varying 
by a factor of a few (see fig0 . This indicates that at later 
times, the information of the initial conditions is completely 
lost. 

More interesting is the way how the cloud evolves. For small 
inclinations (ip ^ 30°), no significant distortions occur and 
the mass accumulation and the final mass are comparable 
to cloud masses formed by purely head-on collisions (within 
factors of 2-3). The inclined flow is aligned with the mag¬ 
netic field very fast. For highly inclined streams the conden¬ 
sation from the WNM to the CNM sets in later due to the 
above mentioned processes. At the same time, mass growth 
is stopped and a short phase of mass loss is evident as a di¬ 
rect consequence of strong shearing motions (see fig. [7]). But 
as soon as the strongest unstable fronts have vanished, the 
cloud turns back to a stabilised state with continuous accre¬ 
tion of matter. Figure [8] shows the resulting molecular cloud 
structure for run B3M0.5I50. Shown is the column density 
along the z-axis, that is, perpendicular to the background 
magnetic field. The black dot resembles a sink particle. The 
global morphology of the cloud is mainly influenced by the 
geometry of the colliding WNM streams with additional im¬ 
pact by the misalignment of the flow. It resembles a sheet¬ 
like shape with trailing arms with the one at the near side 
of the tilted flow being more elongated. This elongation is 
due to the later collision of the flows when the bulk of the 
mass has already been compressed. The resultant motion of 
the cloud yields that the still streaming gas interacts with 
the outer edges of the compressed gas by ’pushing’ it away 
from the actual molecular cloud complex, thereby forming 



Figure 7. Temporal evolution of mass for runs B5M0.5I. ... The 
horizontal black,dash-dotted lines have the usual meaning. 


this observed elongated structure. At this time the flow is 
already too slow to significantly compress the gas, implying 
that the gas in the trailing arm is not able to sufficiently 
cool down by thermal instability. It is therefore not able 
to become gravitationally unstable. These shear flows and 
the resulting occurence of trailing arms are possibly seen in 
observations of e.g. the Taurus molecular cloud, i.e. the non¬ 
star forming low column density arm ( |Alves et al . |2014| >. 

In contrast to the evolution of the total mass of the molec¬ 
ular cloud, the evolution of the stellar (sink particle) mass 
is greatly influenced by inclining one flow. The most obvi¬ 
ous indication is the delay of star formation with increasing 
misalignment (see fig. [9|. Due to the misalignment, mag¬ 
netic pressure and magnetic tension act as opposing agent 
against gravity. In addition, the shear flows disrupt density 
enhancements and thus the transition from the WNM to 
dense, cold structures is hampered. Once, the turbulence 
has fully vanished, the cloud is still subject to its fast bulk 
motion. Clumps within the complex can only grow by ac¬ 
cretion of matter from the immediate environment, because 
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Figure 8. Column density map for run B3M0.5I50. The total 
integration length is 60 pc. The resulting global cloud structure 
is due to the initial compression by the flows. 




Figure 9. Evolution of individual sink particle masses for three 
inclinations of runs B3M.. .1.... The higher the inclination, the 
later star formation begins. In addition, the star formation rate 
and efficiency are strongly influenced (see also fig. [2j. 


of the lack of turbulent compression, which could provide 
the seeds for gravitational unstable cores. The shear due to 
the misalignment also yields a less compact cloud. The ma¬ 
terial is not fully compressed by the two flows. Instead, a 
great amount is at first compressed and enters a phase of 
oscillating motions and dispersion due to shearing motions. 
The accumulation of enough Jeans masses to render the gas 
gravitationally unstable is delayed and also very inefficient, 
since the denser regions are greatly scattered and do not 
possess enough mass. 

As can be seen from fig. [8] there is at least one sink particle, 
indicating ongoing star formation. We stopped the simula¬ 
tion here, because star formation proceeds from there on (as 
can be seen e.g. from run B3M0.4I0). 


4.3 Comparing cloud dynamics in magnetically 
differing environments 

Figure [To] shows the column density in the direction 
perpendicular to the background magnetic field for three 
different initial magnetic field strengths (from left to right: 


\B] — 3 /xG, \B] — 4 /xG, and \B\ — 5 /xG) with an initial 
tilt of ip = 60°. The weakest field case shows a strong 
distortion of magnetic field lines as well as the onset of star 
formation. In comparison, the stronger fields show a more 
ordered magnetic field, which shows no clear deviation 
from its initial uniform alignment. For \B] — 4 /xG one 
can infer some large scale modulation of the field due 
to global dynamics as a resulting imprint of the large 
inclination. The morphology of all three molecular clouds 
is very similar, although some local differences occur. The 
main cloud (having a sheet-like shape) is more compact 
for weaker fields, whereas the difference between the two 
strong magnetisations is negligible. This attribute results 
from the thermally dominated gas. The magnetic field does 
not control the gas dynamics and thus is forced to follow 
the motion of the fluid. Once, local density enhancements 
condense out, the magnetic field is dragged inwards together 
with accreting material. In the cases of more realistic fields, 
it is the magnetic field that dominates the fluid motion and 


that keeps the cloud coherent (see e.g., Hennebelle 2013). 
At the same time the trailing arms now occur to be slightly 
denser. These arms are magnetically supported and thus 
more stable against shear flows and mixing by large scale 
fluid instabilities. 

In contrast to a 3 /xG-field, there is no star formation for 
the cases of 4 juG and 5 juG, yet, although fig. [4] indicates 
that these highly magnetised clouds are also more massive. 
The greater total masses and the lack of star formation 


combine to a picture of a fragmented cloud (see fig. 14). 
Any intrinsically driven turbulence is subalfvenic and the 
magnetic field thus stays coherent. Such a field configura¬ 
tion has also been observed via polarised emission from CO 


(e.g. Li et al. 2010 2014). So, what is the basic impact of 
the magnetic field on the star formation process? As long 
as accretion happens along the magnetic field, the influence 
of the latter can be safely ignored. Once, the gas begins 
to fragment, subcritical regions are produced, as long 
as the parental fragment was only slightly supercritical. 
Thus, accretion along field lines has to continue in order 
to generate supercritical fragments. Otherwise, magnetic 
pressure will drive the gas out of the potential well and the 
fragments stay subcritical and star formation stops. 

The temporal evolution of the line-of-sight component of 
the magnetic field as function of column density is shown 
in fig. mi Different colours or linestyles denote different 
evolutionary stages. The straight lines indicate different 
criticality conditions according to various studies and using 
different approaches for deriving this condition (|McKee 


et al.|p~993l |Shu et aT7||1999| |Crutcher et al.||~2010| [lT et al. 


2014). From left to right the initial magnetic field becomes 


stronger and from top to bottom the inclination increases in 
steps of 30°, starting at cp = 0°. All cases have in common 
that the column density gradually increases at constant 
magnetic field magnitude, indicating gas accumulation 
along the field lines (see also Crutcher et al. 2010). 

The thermally dominated case shows signs of early frag¬ 
mentation and field compression, giving rise to an increase 
of the magnitude at relatively low column densities. These 
effects render the whole cloud magnetically subcritical. 
The clouds then undergo different dynamical phases with 
varying contribution of the magnetic field, i.e. times of 
pure accretion along the field lines, twisting of the field by 
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\B\ = 3/j.G 


\B\ =4nG 


\B\ = 5/j,G 



collapse and compression, and finally amplification of the 
field by large scale collapse and increasing column density. 
In between there exist stages, where the cloud shows signs 
of supercriticality. We here point out that the ordinate only 
shows the average line-of-sight component, i.e. there exist 
indeed supercritical regions that are not significant in terms 
of mass or volume fraction (see 2D histogram). Comparison 
with the weakest field shows that star formation is imme¬ 
diately initiated, when the gas becomes supercritical. The 
collapse proceeds and more material is dragged into the 
potential well. The resulting magnetic field amplification is 
still too low and finally it diffuses out of the central region. 
After the sink particle has formed, some of the field lines 
relax, thereby decreasing the density in some regions. As 
time proceeds the cloud becomes more compressed due to 
its global gravitational collapse. 

For better visualisation, the columns for B = 5 /xG 

and B = 4/xG are shown in the column density range 
19.5 ^ logiV ^ 21 and 1 /xG ^ |B| LOS ^ 25/xG, since 
there occurs no significant amplification of the magnetic 
field during the evolution of the molecular cloud. This is 
indeed very intriguing, because observed magnetic fields 
are far larger in magnitude. We only see motion along 
the field lines, as has already been mentioned before, but 
we do also see no sign of gravitational contraction. Only 
some small modulations are seen, especially in the case of 
the cloud formed by head-on collision and B = 4/xG, but 
this amplification is less than a factor of two and thus not 
significant. At low column densities instead one can infer 
a slight ’global’ amplification of a view percent. This can 
be accounted for accretion of mass from the diffuse halo 
surrounding the dense cloud. These data already indicate 
that the uniform component of the magnetic field is the 
leading component and no clear tangling of the field is 
observed. Furthermore, every process of gas accumulation 
perpendicular to the field lines is instantaneously balanced 
by magnetic forces (see fig. 12). The flow cannot become 
dynamically important in order to bend the field lines and 
to render the magnetic field supercritical. 

Even in the case of high diffusivity or large inclination, 
there is no amplification and/or tangling seen, indicating 
that diffusion processes might play only a minor role in 
rendering the field supercritical. If one takes a look at 
fig. Il3] it is obvious that the gas is highly subcritical. 



Time [Myr] 

Figure 12. Temporal evolution of the ratio of turbulent ram 
pressure to magnetic pressure for three different runs. On average 
the magnetic field dominates and the ram pressure is not sufficient 
to induce motions perpendicular to the field lines. 


Shown are mass histograms as function of the normalised 
mass-to-flux ratio for the three magnetic fields at three 
late evolutionary stages. It is only for the weakest magnetic 
field that the gas shows some sign of evolution. One can 
clearly identify the power-law tail (which can be accounted 
for the iV-PDF) and its growth as more mass enters the 
supercritical regime. In contrast, the higher magnetisation 
cases show roughly no evolution. Once a given distribution 
of the gas has developed it is seen to be globally stationary. 
The difference between the 4 /xG and 5 /xG cases are small, 
i.e. the 4 /xG case develops some larger mass-to-flux ratios. 
However, both regimes are far from being even critical. 

Run B5M0.5I60AD includes the process of ambipolar 
diffusion in addition to an initial tilt. As was already 
mentioned in the remarks of table [T] the simulation was 
stopped at t « 12 Myr. The subsequent evolution of the 
cloud showed no significant difference to the runs without 
ambipolar diffusion. Also in this case, the shear flows tend 
to suppress the formation of dense cores, where ambipolar 
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Figure 11. Evolutionary path of the cloud in N — Blos space. In this case, Blos — B x . From left to right: B = 3)U,G,B = 4p,G, and 
B =s 5/iG, respectively. From top to bottom: <E> = 0°,4> = 30°,<E> = 60°. Different colours denote different times. Note the different 
data range for the stronger fields. Also shown are the criticality condition (Crutcher et al. 2010 Crutcher 2012 dot-dashed black line), 
corrected for projection effects ( |Shu et ah]|1999| solid straight black line), and assuming equipartition of turbulent and magnetic fields 
( |McKee et al. 1993] dashed black line). Colour coded is the mass as function of column density, N , and line-of-sight magnetic field Blos- 


diffusion would be most efficient. 


4.4 Dynamics of dense cores 

Due to turbulence, overdensities occur which become gravi¬ 
tationally bound. Figure [14] shows the evolution of the dens¬ 
est regions within the formed molecular clouds, i.e. of these 
with minimum density of n — 1000 cm -3 . The left image 
shows the evolution of mass, the right panel the evolution of 
thermal and magnetic Jeans numbers. The temporal evolu¬ 
tion is shown for the whole simulation, thus earlier ending 


graphs indicate the complete lack of gas with the respective 
minimum density from this time on. 

For all runs with zero inclination, the compression by the 
two converging streams induces a transition to dense mate¬ 
rial with a few thousand solar masses. But as soon as the 
ram pressure of the confining flows becomes weaker these 
dense regions re-expand, showing that the regions were only 
pressure confined entities. For run B3M0.4I0 a phase of in¬ 
creasing mass follows, which is mainly due to accretion of 
matter along the magnetic field lines. In the end the dens¬ 
est regions of the molecular cloud reach a total mass of a 
few hundred solar masses. Comparison with run B3M0.5I60 
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|S| = 3/j.G |f}| = 4 juG |B| = 5/xG 





Figure 13. Mass as function of the mass-to-flux ratio for B3M0.5I60, B4M0.5I60, and B5M0.5I60 (from left to right). Different colours 
denote different times. The vertical line again indicates the critical value. 



Figure 14. Temporal evolution of the gas with minimum density of n = 1000 cm -3 , resembling the densest parts deeply embedded in 
the cloud. Left: Temporal evolution of mass. Note that for runs B4M0.4I0 and B5M0.5I60 there are stages where no dense material exists 
[indicated by the linear increasing (B4M0.4I0) and the non-varying stages (B5M0.5I60). The former is just the line connecting two data 
points!] between 10 and 20 Myr. Right: Evolution of the number of thermal and magnetic jeans masses. Colours correspond to the left 
figure. 


shows a difference of only a factor of a few in the end of the 
simulation. 

The most striking difference is the first evolutionary phase, 
where the cores of run B3M0.5I60 undergo strong variations, 
because of the additional shearing motions, after they have 
firstly formed at far later times. 

Run B4M0.5I0 already shows the influence of the stronger 
magnetic field. The decrease in external ram pressure by the 
converging flows also induces a re-expansion of the dense 
material within the cloud complex. But now the magnetic 
field is already strong enough to ensure a less efficient mass 
accretion. At around t ~ 10 Myr no dense material exists. 
This stage lasts until t « 20 Myr, where the global collapse 
of the molecular cloud yielded strong enough compression 
to form dense material agairj^] Strong internal variations of 
the cloud then lead to a highly varying mass evolution. In 
the end, masses similar to run B3M0.4I0 are reached, and 
stars start to form. Interestingly it takes roughly 20 Myr 
for stars to form after the reoccurence of dense cores. This 
already indicates that for shearing flows the onset of star 

3 Note the linear increasing interval is simply the connecting line 
of two data points at 10 and 20 Myr. 


formation with magnetic and thermal energies in equiparti- 
tion is further delayed to far later times. But during such 
a long evolution, the clouds would then be subject to large 
scale Galactic processes and our setup would not be appro¬ 
priate. 

Further increase of the initial magnetic field strength shows 
even more dramatic changes in the overall evolution of 
the densest regions within the molecular clouds. In run 
B5M0.5I0 the existence of dense cores ends after t ~ 12 Myr, 
showing the complete lack of unstable cores after the com¬ 
pression by the bulk flows. Although a molecular cloud 
forms, it does not possess any region, which could possibly 
undergo gravitational contraction to form stars. However, 
the first evolutionary stages during the compression of the 
flows shows that the strong fields lead to higher masses of 
the dense gas due to the influence of the field. Inclining one 
WNM stream now shows striking difference. At first, the 
build up of dense cores starts out at later times as in the 
case of run B3M0.5I60. But the diffusive nature of this for¬ 
mation mechanism leads to the build-up of denser regions 
up to the end of the simulation, although there are stages 
where no dense cores exist. The evolutionary track of the 
dense gas is mainly influenced by the cloud motion and the 
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Figure 15. Normalised mass-to-flux ratio as function of 
core mass for runs B3M0.4I0 (circles), B3M0.5I60 (triangles), 
B4M0.4I0 (squares), and B4M0.4I60 (diamonds) at different evo¬ 
lutionary stages. Note that a comparison of equal temporal stages 
is not possible due to the lack of cores with densities of n ^ 
1000 cm -3 . For runs B5M.. .1.. . no such cores were found! The 
mass-to-flux ratio is normalised to fi c = 0.16/VG 
Nakamura 1978|. 


Nakano & 


magnetic forces. The whole dense material is thermally and 
magnetically highly stable, with the latter being the domi¬ 
nant aspect. Thus, although far more diffusive, the magnetic 
field is still able to suppress the formation of unstable cores 
and the subsequent star formation. 

The stability of the dense regions is also indicated by the 
mass-to-flux ratio (see fig. 15). Only if the magnetic field is 
sufficiently weak, gravitational energy dominates over mag¬ 
netic energy and the inner regions of the molecular clouds 
are rendered magnetically supercritical. In case of a 4/xG 
field a great spread of mass-to-flux ratios is observed with 
all being subcritical at the stages shown. During the fur¬ 
ther evolution of the cloud, magnetically supercritical cores 
form. Additionally, by inclining one stream, low-mass cores 
are generated, which have approximately the same mass- 
to-flux ratio as in the case of head-on colliding streams, in¬ 
dicating that the magnetic field diffuses out of the regions. 


4-4-1 Analysis of the densest cores 

Now we analyse the three densest cores in more detail. For 
simplicity, we have assumed that the cores are spherical en¬ 
tities. Figure shows the normalised mass-to-flux ratio, 
the turbulent sonic and the turbulent Alfven Mach number 
as function of radial distance from the centre of mass. 

For run B5M0.5I0 the mass-to-flux ratio is subcritical and 
constant throughout the whole core. This has two implica¬ 
tions: 1) Magnetic support is sufficient to keep the core sta¬ 
ble and 2) there is no evidence for accretion of matter along 
the field lines (which would increase the ratio locally). Be¬ 
sides being higher, the mass-to-flux ratio for run B3M0.4I0 
shows some variation as function of radial distance for all 
three cores. The centre of the densest core (solid line) is seen 


to make a transition to a supercritical state surrounded by 
a subcritical halo. Although the other two cores are subcrit¬ 
ical as a whole, they show the same signature. The cores in 
run B4M0.4I0 show a state between these of runs B3M0.4I0 
and B5M0.5I0. As expected, the mass-to-flux ratio of the 
densest core decreases with increasing radius. The two other 
cores show only a roughly constant ratio. The mass-to-flux 
ratio of the first core is still subcritical, but the transition to 
a supercritical state is achieved at slightly later times. Note 
that the location of the maximum mass-to-flux ratio in this 
core does not coincide with the centre of mass. 

The dynamics of the cores can be analysed by looking at the 
turbulent Mach numbers (see middle and right panel). The 
cores are subsonic and subalfvenic , hence showing that 1) no 
strong compressions within the dense material occur and 2) 
the magnetic field prevents the gas from accumulating into 
denser unstable fragments. This is true for all clouds and can 
be interpreted as an imprint of the initial conditions. Only 
for the densest core in run B3M0.4I0, the outskirts are seen 
to be slighty supersonic and superalfvenic, indicating turbu¬ 
lent accretion onto the core. The difference between the two 
Mach numbers is less for the weakest field, since it is not 
able to fully prevent collapse. If the field strength is higher, 
magnetic tension will accelerate the gas, while relaxing the 
field lines. This is why the sonic Mach number is slightly 
higher, but nevertheless the motions only reach subsonic or 
at most transsonic states. At these densities (n ~ 10 3 cm -3 ), 
the flow seems to be mediated by the magnetic field lines, 
which in every case tends to suppress the build up of turbu¬ 
lent vortices. 


5 SUMMARY & DISCUSSION 


In this study we have presented the results of MHD sim¬ 
ulations of colliding flows with varying initial conditions. 
The strength of the turbulent velocity fluctuations, of the 
background magnetic field as well as the alignment of one 
of the WNM streams with the magnetic field was changed. 
We have shown that dense clouds can form independent 
of the initial conditions, but that their final mass and 
dynamics are mainly controlled by these (see tab. [2|. 
Increasing initial turbulence lead to lower cloud masses due 
to less coherent gas streams. Oblique flows still lead to 
clouds with masses comparable to what has been observed 
recently and stronger magnetic fields will generally lead 
to more massive molecular clouds. The first point seems 


at first a contradiction to Inoue & Inutsuka (2009) who 


stated that for larger inclined flows no dense, molecular 
clouds can form. However, here the cloud accretes mass 
and becomes molecular with time. As can be seen from 
fig- III the onset of the formation of dense gas is delayed 
with increasing inclination. This is indeed consistent with 


Inoue & Inutsuka (2009), because the first few Myr are 
characterised by HI gas with densities below the threshold 
density of n — 100 cm -3 . 

Molecular clouds are able to condense out of the WNM, 
independent of the magnetic field strength. However, only in 
the cases of fairly weak initial magnetic fields, the formation 
of stars could be initiated. Starting with subcritical HI 
flows, the magnetic flux loss is in no cases sufficient to allow 
the build-up of supercritical cloud cores. The tendency of 
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Figure 16. Radially averaged profiles of the three (indicated by different line styles) densest cores within the molecular clouds for runs 
B. . .M.. .10. Left: Mass-to-flux ratio. Middle: Turbulent isothermal Mach number. Right: Turbulent Alfven Mach number. The radial 
distance is evaluated with respect to the centre of mass. For runs B3M0.4I0 and B4M0.4I0, respectively, the data is shown shortly before 
the first star appears. The values for the turbulent Mach numbers mimic the initial conditions, i.e. that the turbulence is still subsonic 
and subalfvenic. Note, the masses of these cores range from w 4 Mq (third massive in run B5M0.5I0) to « 250 M© (most massive in run 
B3M0.4I0). 


Table 2. Typical cloud parameters for gas with n ^ 100 cm 3 


at the end of each simulation with equal turbulent and flow Mach number. 


<t> 

in 

Time 

Cloud Mass 

SF? 

Stellar Mass 

SFE 

E ma g 

Velocity Dispersion 

Njeans,mag 

(°) 

( M G) 

(Myr) 

(10 3 Mq) 

(yes/] 

no) (M©) 

(%) 

(l0 46 erg) 

(km/s) 


0 

3 

21.79 

40.81 

yes 

1346 

3.2 

23.59 

0.99 

23.99 

30 

3 

18.29 

21.84 

yes 

202 

0.9 

12.28 

1.29 

12.51 

50 

3 

18.94 

14.5 

yes 

14.20 

0.1 

6.39 

0.73 

9.86 

50 a 

3 

21.04 

17.9 

yes 

43.83 

0.2 

7.09 

0.84 

14.28 

60 

3 

21.49 

15.2 

yes 

95.68 

0.6 

7.02 

1.45 

10.68 

0 

4 

39.52 

18.33 

yes 

59.32 

0.3 

13.63 

0.53 

5.70 

30 

4 

27.13 

23.93 

no 

— 

— 

21.53 

1.21 

6.30 

60 

4 

27.28 

16.15 

no 

— 

— 

12.29 

1.74 

6.29 

0 

5 

24.64 

43.29 

no 

— 

— 

52.71 

0.46 

6.10 

30 

5 

22.94 

20.39 

no 

— 

— 

22.94 

1.02 

3.07 

40 

5 

37.93 

11.56 

no 

— 

— 

12.08 

1.07 

2.16 

50 

5 

37.03 

17.57 

no 

— 

— 

19.63 

1.42 

3.39 

60 

5 

31.93 

7.31 

no 

— 

— 

6.48 

1.79 

1.55 


a: Different initial random seed for the turbulence. 


the magnetic field to realign itself with the initial direction 
is a crucial factor for the overall evolution. In order to 
circumvent this problem, non-ideal MHD was resembled 
by means of tilted collisions. Increasing inclination leads 
to increased diffusivity of the magnetic field. The variation 
of the inclination as well as the flow dynamics showed no 
tendency for faster accumulation of gas or faster transition 
to thermally dominated regions, since the flow dynamics 
is rather controlled by the appearing shear flows than 
magnetic diffusion. 

We therefore stress the role of magnetic fields in the context 
of molecular cloud and star formation. We point out the 
complete lack of supercritical regions for realistic initial 
field strengths. As was shown in fig. |15[ the normalised 
mass-to-flux ratio ranges from 0.3-1.7 (B = 3/iG) and 
0.3-0.7 (B = 4/xG), respectively. At least the former case 
for weak fields compares well with the results of |Chen| 
&; Ostriker| (20141 /x/Mc - 0-5 - 7.5) as well as with 


observations (r/r c 


Troland & Crutcher 


2008). From 


the observational side, HI clouds may be supercritical as a 
whole, but their observed, dense subregions be subcritical. 
The question remains, how clouds achieve the transition 
from sub- to supercritical. 
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APPENDIX A: ESTIMATE OF THE 
EFFECTIVE MAGNETIC DIFFUSION 

We here give a simple estimate for the dependence of the 
ambipolar diffusion coefficient on the magnetic field strength 
and the inclination of the flow. Specific numerical values are 
not of special interest here. 

The AD diffusion parameter is ?7 ad oc B 2 . The respective 
numerical diffusivity is given by the product of the grid size 
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Time [Myr] 


Figure Al. Shown is the deviation from the induction equa¬ 
tion due to diffusion in the ideal MHD limit for inclinations of 
0° (solid) and 60° (dashed). Different colours indicate different 
initial magnetic field strengths. 


and the Alfven speed. The ratio ryAD/^num estimates the 
influence of magnetic to numerical diffusion. In order to re¬ 
ceive the diffusion for tilted flows, we make use of the fact 
that inclined flows will generate tilted field lines. This tilt¬ 
ing can be interpreted as impact of magnetic diffusion. The 
perturbed field SB (which is here important for the diffu¬ 
sion coefficient) and the initial background field are then 
related by SB = Bosin(ip), where the tilting angle of the 
perturbed field, x, and the inclination of the flows, ip, are 
related by x — 90° — <p. The resultant ambipolar diffusion 
coefficient is then modified to ?7ad oc Bq sin 2 (y?), which gives 
zero diffusion for aligned flows and maximum diffusion for 
the perpendicular case (see fig. Al). 















